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ABSTRACT 


This thesis is a computer simulation of secondary electron 
emission (SEE) from monocrystalline metals and alkali-halides 
using a mod Phen) TaeGHOE of the Harrison, Carlston and Magnuson 
Single collision theory of SEE. Three cases of SEE are investi- 
gated: the angular dependence of SEE from Cu bombarded by Ar, 
the dependence of SEE as a function of energy for rare gas ions 
nomiabiy incident on the (100), (110) and (111) faces of 
Metal Single crystals, and the dependence of SEE as a function 
Saeecnerpgy for Ar and Ne’ ions Nomnally ancident on the (100), 
mieoo sand (lll) faces of KCl. The theory does not accurately 
describe the angular dependence of SEE for monocrystalline Cu 
targets, but does accurately predict the modified sec 6 depen- 
dence round experimentally in polycrystalline studies. For 
ine metal tarcets, the difference between the theoretical 
kKanetic secondary emission result and the experimental datum 


femochtlited as potential secondary emission. The alkali- 


halide SEE simulation agrees reasonably well with experiment. 
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I. INTRODUCTION 


Secondary electron emission is the release of electrons 
from target materials bombarded by primary ions. The first 
observation of secondary electron emission in 1889 was by 
Villard who noticed electrons being ejected from the cathode 
imeeeaischarge tube. Secondary electron mission (SEE) is 
Siamedeterized by the coefficient of ion-electron emission, 
YY, Which is defined as the number of electrons liberated per 
ineident ion, 

Seconaary electron emission can be produced by two 
mechanisms, potential emission (PSE) and kinetic emission 
moby.) Potential emission arises from the fact that as an 
ion approaches a target material it may become energetically 
feasible for a free electron to be emitted from the target 
surface either by Auger neutralization or Auger recombination. 
On the other hand, kinetic emission arises from the actual 
collision between the electrons of the impinging ion and 
the electrons in tite Gaiuceteatron mw sin Pencral. ¥ depends 
primarily upon the PSE mechanism at ion energies below one 
Kev and the KSE begins to contribute at ion energies above 
one Kev. Although this is a fairly good rule of thumb, 
Wetrabtes such as the work function of the target material, 
Tomecharpe, surface state of the target, target temperature, 
and whether the monocrystalline or polycrystalline nature 


Of the target affect the vadue of the KSE threshold. 


10 





Medved et al. [1] have shown that potential emission 

from Ar and Ar° on Mo increases with energy, with the 
PreweokLe at higher energies consisting partly of a KSE 
Gemponent and partly of a PSE component. In general there 
exists a transition region that is not sharply defined and 
Bometuhis reason it is felt that PSE contributes to secondary 
emission in the energy range above one Kev, even though the 
KSE should be the dominant emission mechanism in the one 

to ten Kev energy range. Both types of emission are investi- 


peaued in this thesis. 


JT 








ieee tUDY JOBJECTIVES 


inel964 Harrison, Carlston and Magnuson [2] proposed 
feeunegle collision theory of secondary electron emission 
from monocrystalline target materials in the one to ten Kev 
Peimary ion energy range. The original paper develops a 
theory applicable to monocrystalline targets bombarded by 
[mierc ¢a4S 10n beams directed normal to the (100), (110) 
mee bll) crystal faces. This thesis will extend this 
theory to explain the following phenomena in the one to ten 
Kev primary ion energy range: 
Mee occondary electron emission from a single crystal 
plane rotated arbitrarily about a crystallographic 
AALS i 


2. Secondary electron emission from metal single 
crystals, and 


3. Secondary electron emission from alkali-halide 
monocrystals. 
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Pit COECONDARY ELECTRON EMISSION THEORY 


Secondary clectron €jyection mechanisms are either of 
meempotential or kinetic type. Potential ejection can be 
caused by resonance neutralization, resonance ionization, 
meer neutralization, or Auger de-excitation. Kinetic 
Pyleectilon requires the direct transfer of the ion's kinetic 
emerey to the electrons of the target atom. Both mechanisms 
are discussed below. 

Resolance Neutralization occurs when it is energetically 
measible for an electron in the target material to "tunnel" 
Bpipeugh its potential barrier and approach the position of 
miewricoming don. ihis stage of the process increases the 
energy of the ion and target atom system by an amount equal 
boethe energy needed to remove an electron from the metal. 
This energy can range from near zero to $9, the work function 
of the target material. In the second stage of resonance 
iewiralization, the electron and the ion combine, effectively 
mesiniLinity, to reduce the system energy to an amouht Vis 
Pies ronization potential of the ion. This mechanism is not 
g@urte aS simple as it seems, eigce lt iswat.ected by the 
amace potential and coulombic forces of the electronic 
PenUctNne OL the atoms in the crystal. According to the 
FranckCondon pHmieiple. the electronic transition of the 
ellecrrometrom the target material to the ion can only occur 


with appreciable probability when the position and velocity 


eS 





ieee, 10n are not affected by the transitiem. By this 
Mmeaiciple we expect a crossover of the initial and final 
potential energy curves of the system. Taking the image 
meecential into account, resonance neutralization can occur 


for reasonable separations of ion and target atom when 


(V. - ie) - o> + ile < Ws - o, 


and 


(Ve 


go Me) a 


Therefore do < (V. - Mead 12 te < ie 


> 


meme ¢ 1S the work function of the target material, V;> the 
Homization potential of the ion, ee tne excitation potential 
of the primary ion, I; the surface image potential and Wo> 
Baemitttial Kinetic energy of the primary ion. 

imessecond possible potential ejection mechanism is 
Mmesolance 1Onization which is essentially the reverse process 
GiewesoOnance meutralization. In resonance ionization an 
--erecapatom at Infinity 1S lonized into an electron and an 
won. - the electron then proceeds to the metal surface and 
@ceupies any Of the unfilled available states above the Fermi 
level of the metal. As in resonance neutralization, there is 
a range of ion-target atom separation for which the process 
immostelikely £0 occur, There iis also the possibility that 
a highly excited atom will be ionized very close to the 
target surface with a subsequent high probability that a 


Secondary electron will be produced. 
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Pithoweh the above processes are possible, they are 
Moeretie ones that make the largest contributions to secondary 
electron emission at low energies. Auger processes, which 
Secur at distances of approach of 3 R or less, are the 
primary mechanisms contributing to PSE in the low Kev 


energy range. These can be summarized as follows [3]: 


Puppet aneWwtira lization: 


eens ne, >X+e + (n= 2jie" ; 


and 
Nitec e-EXCl tation. 


m - ~- _ 
ee nee X + e@ + (n-lje, ; 


where X and X" refer to the bombarding ion or metastable 
atem, respectively, and ne, iomtne mumver OL etree ellectrons 
Mititaniy sin the target material. In both cases, the final 
Seate Of the system is such that a free electron and a 
mMeutral atom are produced. Auger neutralization occurs 
when the energy of the initial system, target and ion, is 
greater than the resulting system, X +e + ive ie. 5 Bie 
infinite separation. Schematically the process is shown in 
Preure Ie 5]. Auger de-excitation, on the other hand, 
involves a metastable atom which approaches the target 
Material, As a result, the initial energy curve, shown in 
Pecinewz|Siaes mot modified by an image potential and 


DemGcCeetnesiiitial and final.energy States never cross. 


is 





Eke nay Here +n-2)e ,, 


FE NERGY (Kev) 





| 2 Se aaerErs - 
S(A) 


Figure 1. Potential energy diagram for a helium ion as it 
undergoes Auger neutralization at a tungsten 
Ssurrace (3). 
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ENERGY (ev) 


be 7 Ne+e +(nDe 
40 | e (rn ey 





PeirCmzaet Otential energy diagram for a metastable atom (He) 
as it undergoes Auger de-excitation at a tungsten 
Surface (3). 


ds 





Mere itwO processes, tO a preat extent, account for low 
energy (ev range) secondary electron emission. 
werkewenrssion 15 the result of the direct transfer 
of the incoming ion's kinetic energy to the atom's electrons. 
femteor this energy is absorbed in the recoil of the aton, 
while some of the remaining energy is transferred to the 
emeetrons Of both the ion and target atom, causing an 
Peemcation Of both of their electronic structures. As the 
meemm@uclei: Separate, some of the perturbed electrons are 
Imeievated aS secondary electrons. The ion and recoiling 
Meet atom are slowed down within the crystal structure 
Gertie target material and are likely to undergo Auger 
Mmemetralization or de-excitation. If this process occurs 
weertemently close to the surface of the target material, 
m@iere iS a probability that these electrons also will 
eceane from the target materials surface and contribute to 


kinetic secondary electron emission, 
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iv) THEORIES CONCERNING SECONDARY ELECTRON EMISSION 


Pee tHEORY OF PARILIS AND KISHINEVSKII 

Pemalreetie theories proposed to explain kinetic secondary 
Peectron emission, the theory of Parilis and Kishinevskii 
(PK) seemed to be one of the most promising [4]. PK 
fevided the secondary electron emission problem into two 
separate subproblems which were analyzed separately and 
Bee cquently collated to develop the theory of KSE. These 


Exo subproblems included: 


ieee ine search for a sufficiently effective mechanism 
SeelreanotersoL tne kinetic e€nergy of a moving 
particle to the (target) electron; and 


Pine aescription of the yield of excited electrons 
from the target. [4] 


Mieassumed that the solution to the first problem was of the 
femme proposed by 0. B. Firsov for inelastic scattering of 
Hon-atom collisons in the Kev range [5]. Secondly, PK 
faeWalized tie secondary electron emission process as one 
toewhich the primary ioimelomlses Manyoatoms as it progresses 
pillonp sits trajectory in the crystal. Electrons produced 
in this manner, Ob Veiurther wonization of the primary 
fone then proceed to the Surface ot POE collisions and 
energy llosses as welll as possible capture in route. The 
combination of the electron production phase and the elec- 
tron escape phase then yields the KSE from the crystal. 

The actual collision between the ion and the atom was 


seen by PK as being accompanied by an overlap of the 


Is, 





emecetron shell structure of the two collision partners, 
Guring which time there occurred an electronic interaction. 
Pierey was transferred between colliding particles, exciting 
their electrons to a higher energy state. To calculate this 
meercation energy, a statistical picture of the electron 
distribution of the particles was used. The theoretical 
relationship, based on the Thomas-Fermi (TF) model of the 
meom, leads to the following formula for the energy trans- 


f/rrea between the colliding particles: 


mo? 


ann? 





tr 
I 


f (f, 9° ds) R+ aR, 


iene R is the radius vector connecting the centers of the 
Sgmliding particles, S 1s the surface separating the domain 
Giesetion tor the potential of the two particles, and 4 
Hoeeve potential on this surface given by: 


ys 
erie Ih ILS (7h ie) as 


¢= [4-1 x | . es 


Me 


where Z. isetie charge on Gach particle, arp = h?/me? and 
faceetnen lr seneening function [8]. Physically, the inter- 
fetioumcan be thought of as analogous to the collision of 
Pvomtenniseballs. As they impact, both deform, with the 
vector joining their centers being perpendicular to the 
area S of the plane of interaction which separates them. 
Obviously the total energy transferred is a function of the 


relative velocity with which the ion and atom impact and 
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the electron distribution associated with each. In terms 
@recne tennis ball analogy, the greater the relative impact 
Hetocity, the larger the deformed area S, the shorter the 
Geaiis vector distance R, and the greater the number of 
efectrons involved. 

The analysis was rewritten by PK in terms of the impact 
Bawameter associated with the collision. Hence by changing 


Variables 


Gem eels = VAT) 


Zz 
V(r) _ P 
gales tie a7 


E(p) = —5 (2 + Z dR 


ye X{{p) do 

R/2 Ne) 
where Uy im—tneominitraleyelo@crty of the ion, E is the 
Peeitation energy of the colliding particles, p is the impact 
parameter ander) as) the potential of repulsion between the 
particles. 

MOmObtainethne correct behavior of the emission curve, 
PinemGross section of electron ejection into the conduction 
band had to be calculated. Because more than one electron 
iayebeslitted trom the valence band to the conduction band, 


the cross section was determined to be 


Care enCeonstanmt «(Probability of ionization) 


P -(Impact parameter) 


Zi 





which in the integral limit approaches 

Je 

C= 27 f i 280 dp 

O 
Here Py is the impact parameter at which E(P,) = 6 - 9, 
fiere 6 15 the energy depth of the filled conduction band, 
fetus ctie work function of the target material and J is the 
average ionization potential for the outer shells of the 
m@eom. For 1/4 < Z,/Z. <4 PK found the cross section to 


be given by 


2 
ieee (eae) 


OU ——— a as Ys ana W se ee Oe, 
& es) 


and S(uo) is approximated by 


7 


S(u,) = 5.25 u, arctan (.6 x 0 Gis ce Ts Aa ee 


O min 


where u isechnestiresiold velocity. “After determining 


min 
Ene probability that an electron will be promoted from a 
filled band to the conduction band, PK assumed that most 
Sem@iesceclectrons remained in the conduction band because 
of the surface potential of the metal target. However, 
Peeuisewaenole vis left in the valence band, an Auger recon- 
bination iS possible with subsequent emission of an electron 


Preoneciemmetalwes Ihe probability that this process occurs 1s 


given by PK as 


w(6) = 0.016 (8 - 26) 


ZZ 





PK assumed that both PSE and KSE were due to Auger 
Pocesses [4]. 
miter the Secondary electrons have been produced, they 
fweerer CcOllisions as they proceed to the surface of the 
metal. These collisions degrade their energy and increase 
the probability that they will be lost through recombination. 
PK assumed that the number of electrons which reach the 
SUrface is given by an exponential law of the form exp (-x/A) 
Milerye x 1s the depth at which the electrons are formed in 
f@ijemenystal and 4 is the mean free path of the electrons 
Mimete Material. Specifically, the number of electrons 
Eile ted per ion incident on the surface is given by 
y = x (Number of atoms Ce Ober end eanec LOCEron 
x pemacubic. centi - Will be ejected into 
meter) the conduction band 


and a hole left ain 
valence band) 


W(Goopwmbiac. an electron 
will escape to the 
Suigrdee vou the crystal) 


Vinten in the limit poes to 


X 
aC eNeo Gai wi(ohwexp (-x/\) dx , 
O 


where Se lowene maximum depth at which the ion retains 
sufficient energy to ionize a target atom. Upon taking 
into account the degradation of energy of the primary ion 


Foettesuirers multiple collisions and the diffusion cross 


ZS 





Beaton Of the eléctrons, PK arrived at the final expression 


mor Y: 
fe cto}, [or ue ji: Ae (ud) : 
where 
Ao (u,) = exp (-uo/ (KA) fio exp (u2/ (KA)) a) du 
and 
C2 t8m Ns agp 2 ee 
(M, + M,)(z,1/? + 20/2)" ; 


ies: kK theory predicts the following trends [3]: 


Dee Low velocity region: uu. = u 


5 eee Ao (u,) term is 


important..., y increases slowly with Uy» 


(o (u,) - Ac (u,)) = lel =pallo, 2 Ul ye. Hence 


O 
initially yous, ice y 1S ad lanear function 


min 


of energy... 
2. High velocity region: (o(u,) - Ao(u,)) = 


u tan”? (.6x10° /(uo-u eins express jon 


min 
asymptotically approaches tLieestrarcnt. lune 
(GbE e(u Rou) and therefore increases 
Mineamiyewitiie velocity, Extrapolating the 
lineaG, portion back to o(u,)=0 devalue Of the 
IME We Oe u,=1.05x10/ cm/sec is obtained; 

Pio mEsmIndependent Of 10n target combination. 


Meveny site velocities: Here ion penetration depths 


are greater so that electrons are formed further 
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from the surface. The number escaping to the 
surface and therefore the yield is expected 
to pass through a maximum. This is confirmed 
experimentally. 
4. No yield dependence on ion charge is predicted. 
wee aGependence of yield on the ion target combination 


L/2 Ay, 4 


Ot the form ((2,+2,)/ (2, +Z, nic for heavy 


ions and (C2, 
ie the term Ao (u,) Pop OLE Ants eackie \ helde 1S 
independent of ion-target combination to a first 
approximation. 

Gemebependence of yield on the angle of incidence 
Wms eexpected to be of the form Vet eree BD since 
Phempnobdhality Of electron eScape 1S a function of 
MicmslGimecscteaistance to the surface and the proba- 
peRtty sor tOrmMatyon Of electrons is a function of 
Dicmactual distance traversed. 

jemenewetfect of the 1Sotopic mass on yield is 


Mitempated to be connected with the different 


hetardatwonerates tor the isotopic pairs." [3] 


MihTOr stnese conclusions are approximately confirmed by 


the experimental results. 


B. THE THEORY OF HARRISON, CARLSTON AND MAGNUSON 
The PK theory uses many simplifying assumptions which 


introduce approximations into the results obtained. In 
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Particular, the PK theory predicts secondary electrons of 
Oily average energy, two to five ev, and data obtained 

by Wolff [6] of SEE from copper bombarded by ten Kev argon 
ions indicates some secondary electrons emitted have energies 
Preater than 140'ev. 

In the PK theory secondary electron emission was viewed 
eemaecwO-Step process in which ions were first produced 
along a path in the target material and subsequently diffused 
Womene Surface with possible scattering or Auger emission 
weeurmnime along the path. This model does not yield the 
Mmeen-energy secondary electrons which are reported experi- 
femecaiiy. In an effort to explain this discrepancy between 
Eeterimental and theoretical result, Harrison, Carlston 
Boomvaenuson |2| (HCM) proposed a "single collision” theory 
of secondary emission. This theory, though in many respects 
camilar to the PK theory, has as its basic philosophy an 
electron production mechanism which differs from the model 
proposed above. 

HCM proposed a model in which an ion approaching a target 
SierdcestIrst Undergoes resonance neutralization and subse- 
quent ly eomcdes Wit ticmbirstelayer of atoms im the target 
Naterial im the Un-1onized state. Even though the primary 
Particle 1s netitralized before the collision, the terminology 
Smelom etserctained here to distinguish between the primary 
particle and the target atom after collision. For simplicity 
H@Meassumed that the electrons liberated due to the inelastic 


energy transferred between the electronic structures of the 
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ion and the atom all come from the ion, which allowed them 
memptonore the energy partition between collision partners. 
Theoretical work by Harrower [7] indicates that the Auger 
Beoeess requires a time on the order of io seconds while 
mimtenhe one to ten Kev energy range the collision process 


eee mune mini ica es 


fakes On the order of 5 x 10 
Wie statistically, the ion usually does not have time to 
Mewrtralize before its next collision and for all practical 
biaewoses Cannot contribute further to the secondary electron 


emission process. These results lend credence to the single 


Comision model. Mathematically HCM proposed: 


(hk) 
(hk2) S 
Y eM Monee (sna) PS (5) ds", 
KSE 0 
(hk2) 
where y fomencenumoen Of Clectrons emitted per incident 
KSE 


ion when the ions are normally incident upon the (hk&) surface; 
Ss is the impact parameter of the ion on the atom; Err Ls 


ener inelastic energy transfer between the ion and the atom; 


(hK2) 
ne (S,Eqrp 


produce ; ane Eee isS)the probability that the impact 


) is the number of electrons this energy will 


Batameter s Will occur in the (hk&) surface. This equation 
contains three adjustable parameters, the constant K and two 
Sonstents which enter into the calculation of ¢, the electron 
density, and the interaction potential V(r). By symmetry, 
these last two constants must be the same for all orientations 


Ofettespanticular crystal plane, (hk2), in question. 
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Peeording to HCM, the constant K is a factor which accounts 
for three physical processes: 


ie the probability that an electron will be emitted 
Meo the surtace (a value close to 1/2); 


Ze the probability of additional emission by an 
Auger process; and 


Sealine probability that a high energy electron 
escapes before it iS scattered back into a 
Sonduction band. [2 | 
feeleor these processes are orientation independent. The 
factors which appear in the formulation of the HCM secondary 
electron emission equation are discussed in detail below. 
imine Inelastic Energy Transfer 


Inethe HCM theory E Teette ine las tile seller oy 


TEP 
Peemererred in the collision process and is basically the 

Soamewas E(p) in the PK theory. However, as HCM point out, 
because the interaction involves atoms spaced regularly in 


amlttice array, the inelastic energy integral does not have 


@neintinite Upper limit as was assumed by PK. Specifically, 


Umax 
(1-V(r)/E_) dr 
ErpplSsE,) = Kopp - I 
Doe 
(1-V(r)/E,-s?/r*) 
D tae 
min 
TEM 
‘ 9? touade 
5 
2 
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where 


h 
Kopp 7 =e (Eg/ (2m) 42 (25425)7 
Oa 


fee hi-ane-) = 529 2 


H 


and Eo Hestie energy of the ton in the center of mass system: 
zy and Z5 are the ion and atom nuclear charges, respectively; 
eS ms the distance of closest approach; and » is the 


uiemes-Fermi (TF) screening function: 


ec o4 aos 


o(o) = [1 + (a,p) ] 


where 


1/3 
_ App ((2y+2,)/144) 
F° 0.8853. 


and arr 1s one of the adjustable constants introduced above. 
Pyeexperimenting with different sets of potentials, in 
Conjunction with the screened electron distribution, HCM 


found that the best agreement with experiment was obtained 


when the potential V(r) was of the Born-Mayer form [2]: 
Wie) = exp CAtBr) = AY exp (Br) 


with B always negative, and where A is the second of the 
above mentioned arbitrary constants. HCM referred to the 
mibccralecontainane dp as the electron integral. The con- 


vention is retained here. 
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Several important results were noted by HCM in the 
noamilation of their original model. First, the dynamical 
mpepral was extremely sensitive to its upper limit eee , 
whereas the electron integral was not sensitive to varia- 
tlons in Tou: For this reason Tey Was Set at a value some- 
what larger than the nearest neighbor distance Be in the 
Mievotal lattice, and retained this value for all orienta- 
Inmemis OL the crystal lattice. On the other hand, the upper 
omic Of the dynamical integral was extremely sensitive to 
the nearest neighbor distance. Such a behavior is expected 
because as each ion collision with a particular atom termi- 
nates, another must begin with one of its neighbors. From 
a particle dynamics point of view, HCM found that by 
treating the actual collision with the lattice as a binary 
Gemilision between the incoming ion and any single atom 
Pomein the first repeat distance of the crystal plane, the 
theory yielded results which approximated those obtained by 
experiment. Figure 3 indicates the geometrical relation- 
ship between temetayhe quantities associated with the 
G@yaamical integral [2]. 

Pea co litsion) Parameters 

Beeawse @: the repetitive mature of the lattice, or 
iaterce plane, 1t 1s not unreasonable to think of a "repre- 
femeatlve area swhich Characterizes the lattice plane as a 
wiele and hence ats distribution of impact parameters. If 
no lattice were present, and only one target atom, there 


would be no characteristic area and the occurrence of a 
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collision with any impact parameter would be possible, 
fuerding an infinite collision cross section. As an ion 
approaches a plane of atoms in a crystal, however, 

Emere CXists an allowable set of impact parameters for the 
Sieeeon collision. Specifically, in a lattice there exists 
an impact probability distribution, eee Oi eae 
Siystal Orientation. For normal incidence HCM proposed the 
@veas Shown in Figure 4. These areas are representative of 
Elie formal iipact Ci an 10N with any atom on the specified 
Uryoctal surface, viz., by symmetry operations they fill plane 
space and are representative under any of these symmetry 
Smeerations of what an ion "sees" as it approaches the planar 
mimeiace normally. There are many distributions which could 
bemeonsidered probability distributions. HCM considered the 
following: 


meee ihe distribution of values of s measured to the 
Mearest SUurtace atom; 


wei iesaistribution of all values measured to all 
Scousewien project their area into the repre- 
sentative area; 


Pee incedistribution of smallest impact parameters 
wien all atoms are considered [2]. 


Wie third case was the only one which yielded results in 
agreement with the experimental data. The distributions 
for normal incidence on FCC and BCC metal targets are shown 
in rreures 5 and 6. 

Operationally, the distributions are determined as 
follows: First, the representative area is chosen. 


"Points" are shot at the representative area, and the distances 


yA 
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Representative areas for normal incidence. 


Reet (2). 


Figure 4. 





NORMALIZED PROBABILITY DENSITY 
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Distribution: 
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igure 5. Distribution of impact parameters for normal 
° incidence on face-centere i 


d= ciubaemervystals. 
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Between the impact point and the centers of the atoms lying 
in the representative area are measured. For each atom the 
@eevances for points that are measured to be closest to it 
Mu@emtallied. After a sufficient number of points are tallied 
to yield good statistics, the ree ninien for each atom is 
Memented by a factor to account for the fractional portion 
Oreatom lying in the representative area. The atoms! dis- 
Meueucions are then added and scaled so that the total 
wecaeunder the resulting curve is unity. 
3. Number of Electrons Emitted--Russek Model 

One of the major shortcomings of the PK theory was 
that it used average ionization energies to obtain an ioni- 
Zoom cross section. This is to say the PK theory did not 
eomemder the electronic structure of the target atom. HCM 
adopted a model proposed by Russek [9] to obtain the number 
eee 


Semeobectrons liberated per collision, In 


Epp): 
faemOriginal paper, Russek proposed two ionization models of 
mimeerest, the uniform ionization model (UJ), and the staggered 
ionization model (ci. The UI model assumed that as electrons 
were liberated by the collision process, the amount of energy 
required to liberate any one electron was the same as that 
Heauiredsto liberate any other electron. On the other hand, 
the SI model assumed that each electron liberated required 
more energy than was required to liberate the preceding 
electron. For the SI case Russek shows that the probability 


that n electrons will be liberated from a shell, which 


Gontains M electrons is given by: 
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ee CL) ts a 


(192 Gin n Dele 
n 1=0 i 


ION _ , ION ,, ION 


where k < (E E 
=. | TER! en n+l ntl 


ike land Mens the munpen 


Ofeelectrons available in the outermost shell of the atom; 


Eo oN Memrle energy required to remove the Ae Gi these 
n 
electrons and E is the energy available to cause this 


TFF 


memoval, Iwo things should be noted. First, the impact 
Mmeaeeamnecer enters this calculation through the inelastic 


ememey transfer. Second, when Boo = ys one obtains 


nee IE n 


tie Ul model. After normalizing the function, Bes Vile 


z - p (M) = ] 
n=0 n 


then, from probability arguments, the number of electrons 
liberated is: 
oe ee! (M) 
n, (Eppp) = - a n 


n=0 e 


HCM were unable to obtain good results with the Russek 
Mienodely In all cases, y increased much too rapidly as a 
function of energy. On the other hand, the SI model yielded 
ION 


reasonably good results, where spectroscopic values for En 


were used, 
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4. HCM Results 

HCM investigated the Ar - Cu system in detail. 
They first examined the BORE dynamical integral for various 
iypes Of potential functions and found that the PK potential 
which was of the Thomas-Fermi-Firsov type was too hard for 
small values of Eo? Viz., the hard core radius of the target 
Reon Was too large, in this energy range, to yield accepta- 
ble results. Attempts at softening the potential by decreas- 
ing aon failéd. For this reason HCM assumed a potential 
ivietrton of the Born-Mayer type. The potential function 
used in each case was found by matching the exponential to 
the TF screening function at some separation. 

HCM assumed that the experimental data contained both 
CONETIbDUGIONS to the 


potential, and kinetic, 


YPSE? YKSE? 
G@ueaieesecondary electron emission coefficient. Therefore, 
Metieidentified the y at which KSE begins to contribute as 
Ypsg and assumed this value remains constant for all energies. 
iitse value for Ypsp Was then subtracted from the experi- 
Memealiy determined y for the (111) surface at 10 Kev. The 
Gil) value of YKSE at 10 Kev, obtained from theory, was 
Saeed to this point, and all other points were multiplied 
Peeene scale factor used in performing this operation. 
Finally the constant value of Ypsp Was dddedetorali the 
Semicdeticoretical values, The result is the HCM value for 
the total secondary electron emission. The results of the 


Hebieci@eory for the Ar’ - Cu evystcimare presented in Figure 7. 
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Ar* - Cu(SINGLE CRYSTAL) 
NORMAL INCIDENCE 


© (111) FACE 
A (100) FACE 


[) (110) FACE 


1 Z 3 4 8) 6 7 8 9 10 
ION BEAM ENERGY (Kev) 


Ar - Cu results obtained from the HCM Original 
model. Taken from (2). 
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In general, HCM's theory agrees reasonably well with 


feperimental results for all systems studied except Ar - Mo. 
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Vee VODITFICATIONS TO THE HCM THEORETICAL MODEL 


fee LNTERATOMIC POTENTIAL 

Independent of the philosophy behind the basic interaction, 
mene HCM theory, like the PK theory, is highly dependent on two 
Seectric assumptions. First, both assume that the inelastic 
Sl@erey transferred 1s of the form proposed by Firsov. Second, 
both theories must assume an interatomic potential function 
which is consistent with the electron distributions of the 
Homeand atom. Since the Firsov model is presently the only 
Geamerehensive model available for the inelastic energy 
transfer, it has been retained here as an approximate model 
Of the inelastic collision process and is assumed to be 
fieemiciently accurate for our purposes. However, the inter- 
Memron potential used in the Firsov model need not be of 
the Thomas-Fermi form, but may be derived from a more 
ieetrate representation of the electron distributions of the 
ion and atom. 

The interaction potential is derived from the electron 
Gistrtbutions of PRL cole sian Partners. imese potentials 
exmsteon three general levels of sophistication: (1) The 
Bohr and Born-Mayer potentials which aap derived from 
pilce@eron distributions which provide screening of the 
mucieus in a self-consistent, but not necessarily accurate, 
manner; (2) Thomas-Fermi and Thomas-Fermi-Dirac potentials 


which are derived from statistical pictures of the electron 
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Peetr butions of the collision partners; and 
(3) Potentials generated from electron densities derived 
mamerically from Schroedinger's equation. 

Parilis and Kishinevskii used a Thomas-Fermi potential 
mmecion in the:theory they proposed. Statistical models 
Seumece atomic electron distribution, such as in the TF 
model, are based on studies of inhomogeneous gases and 
approximate solutions of Schroedinger's equation consistent 
feeemeelis picture [10, 11]. They require that the total 
electron energy be a minimum and that the electron density 
/mieeerate to the total number of electrons in the system, 
Approximations may be made to evaluate the exchange and 
correlation energies. When the exchange and correlation 
energies are omitted and the kinetic energy is approximated 
Dae plane Wave, the TF atomic model results [11]. If the 
exchange energy is retained and the kinetic energy 1s 
approximated by a plane wave, the TFD atomic model results 
[11]. The TF model of the atom has been used extensively 
by Firsov, whereas the TFD model appears in papers by 
Abranmson [13] . | 

Potential functions based on Hartree-Fock solutions of 
the atomic Schroedinger equation have been discussed by 
Wedepoh! [12], Harrison [13], and Wilson [14]. These 
potentials assume the adiabatic approximation, use 
spherically symmetric electron distributions which have 


shell structure, and take into account the correlation, 
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Fame tic and exchange energies through various approximations. 
M@itey appear to lead to improved calculations of the interac- 
Gmom potential [8, 14, 15]. The present investigation 
memiizes an interaction energy calculation similar to that 
proposed by ved Swaan (2 eeeedoOWweVer. the electron densities 
fPeom which the potential is derived are calculated from the 
Hartree-Fock-Slater self-consistent wave functions calcu- 
meced trom the Herman-Skillman [16] computer program. The 
Header is referred to the original book [16] for detailed 
imrormation on the calculation of the atomic electron dis- 
Mmermpution functions and to the work by Torrens [8] for a 


meneral discussion of the interaction potential. 


See tHE BASIC INTERACTION 
Herrison, Carlston and Magnuson assumed that all of the 
eeeoudary electrons produced in the single collision model 


came from the ion. This assumption was tested by choosing 


various values of Been in the Russek model and noting the 
@ependence of n,(Erpp) on ee This experiment was car- 


+ 
mecdmout for the Ar - Cu system. The results were most 


tevoerable for ao 


equa letomlons evaernent rst lomization 
Hetention of argon, and least favorable for neutral copper's 
fifet 1onization potential. This approximation was an 
attempt to simplify the model so that detailed models of 
energy sharing mechanisms could be avoided. In this 

present investigation, this assumption has been relaxed to 
include contributions to the secondary emission by the 


Eaneet atom. 
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Theoretically, the interaction is viewed as the collision 
of two atoms having spherically-symmetric, shell structured 
Seeetron distributions. As the atoms collide, they 'are 
assumed to form a "molecular" structure whose electron 
momizdtion potential distribution is comprised of the 
SUdtate 10nization potentials of the electrons of the atoms 
mimolved in the collision. Since electrons are ejected on a 


Seer energy’ principle, the values of a 


used in the 
Russek model have been modified by include both primary ion 
and target atom ionization energies which are arranged in 
miereasing order of magnitude, irrespective of the atom to 
mien they belong. The mean number of electrons is a compli- 


Gated function of the inelastic energy transfer and the 


Pieray necessary to effect each electron's removal. 


eee wie DISTRIBUTION OF IMPACT PARAMETERS 

To accomplish the objectives of this present investigation, 
Omewadditional modification to HCM's original paper had to be 
made. Since the angular dependence of y was to be studied, 
a new distribution of impact parameters had to be calculated. 
Fortunately, according to the HCM paper, this modification 
to the mathematical model is the only change required to 
take into account rotations of the crystal plane around an 
Paperary crystallographic direction. No other factor 
entering into the secondary electron emission equation is 


Gerentation dependent [2]. 
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imoepossible solutions to this problem are proposed 
Metow. First, a crystal plane normal to the direction of 
mierdence Of the i0n beam can be found, from which a repre- 
Peleative area can be constructed. Such a formulation would 
memanalogous to cutting the crystal in such a manner that 
imempresented a planar area normal to the incident ion beam 
and finding an area on this plane which would be representa- 
memewor that orientation. Alternatively, a representative 
area constructed from the crystal face which was experi- 
mentally rotated can be used and mathematically rotated 
a@peut the specified crystallographic axis. In practice, the 
meese method 1s computationally restrictive because it 
requires that a representative area on a plane normal to the 
fierdent 10n beam be found for each orientation. For this 
meesonrit iS not pursued further here. 

The second method mentioned above yields a result 
Gideeis More in line with the reasoning behind HCM's 
Ommeinal definition of the distribution of impact parameters. 
Practically speaking, a representative area is found for 
etemecrystal face which fills plane space and is representa- 
tive of what an ion sees as it approaches the crystal plane 
fuoneanearbitrary direction. Unlike the representative 
Mmeicoecor normal inciden.e, angularly dependent representa- 
tive areas must be truly representative for all possible 
rotations of the plane and not simply for one particular 
orientation. The representative area for the FCC (100) 


surface which can be rotated in either a positive or negative 
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sense about the <110> crystallographic axis is shown 
merreure 8. Note that atoms which lie within the first 
mepeat distance of the crystal plane are eligible target 
meom> in the representative area by HCM's definition of 
representative area. Also note that these atoms from the 
lower atomic plane appear to move as the crystal plane is 
rotated in a positive sense about the <110> axis (see 
pueure 9). This behavior can be verified by physically 
metating a stick and ball model of the lattice structure 
about the indicated axis. 

Prec s tne representative area is chosen, its rotational 
properties and corresponding distribution of impact 
marameters are calculated as follows. The representative 
faeomts placed in the first quadrant of a rectangular 
Sooreinate system (x,y) (see Figure 10). In this system, 
the crystallographic axis of rotation is specified by two 
parameters; 

imine angle of the crystallographic axis makes with 
the x-axis of the rectangular coordinate system 
(bey) and 


feeiiie y-axis intercept of the crystallographic axis 
iMierneunectanaular coordinate system (x,y). 


This system of coordinates is then transformed into a primed 
meetanoiiar reterence frame (x',y'} for which the crystal- 
lographic axis of rotation is the x'-axis, and whose origin 
iewtivespoint at which the crystallographic axis intercepts 
the y-axis, (Xp>Yp) (see Figure 10). The equations relating 


milemtitanstOrmation of points, from the (x,y) coordinate 
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Figure 8. Representative area of a (100) FCC crystal face 
mOerwecdeapout <Ul0> im 6¢ither a positive or 
negative sense. 
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Figure 9. Representative area of the (100) face of an FCC 
GCivotaieroOtated im a positive Sense about the 
Sihves darection. 
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Poeureml On Rotation geometry for a (100) face of an FCC 
Crystal about an arbitrary crystallographic 
direc ton. 
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System to the (x',y') coordinate system are: 


<= Kk COS (8 = o>) i (ypsin o + peo o), 
voe= R sin (8 - o) + (x,Sin OS RS >). 
a R = (x? i‘ ee 


6 arctan (y/x) 


pide? is the angle that the crystallographic axis makes 

With the, x-axis. If the system is rotated through an angle 
a about the WeooxiswediStanees to any point (x';y*) appear 
bembe shortened by an amount cos a in the y' dimension, 
leaving the x' dimension unchanged. Also, if the atom whose 
Pmonuected poSition lies at (x',y') before rotation is an 
atom whose actual position is Ro units below the plane being 
momotea, then its new apparent position on a plane perpen- 
Siewlar to the direction of incidence after rotation is 
Mmiveneby the doordinates (x',y' cos a - R,Sin Oe 


Peneral then the final transformation equations become: 


Mame Reeos (6 - ¢) = (ypsin do + Xpcos ta). 
ceeeee Sin (0 - >) + (xpsin do - Y pcos o)) cos a - R sin Q, 


iiewseometry implicit in these calculations is shown in 
Preure 10. : 

iitemaetualecalculation of the distribution of impact 
parameters is performed as follows. A mathematically con- 
Serieted, rectangular grid is placed over the area and 


pommtes which lie within the representative area are 
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systematically chosen. Using the transformation equations 
and the known position of the atoms in the representative 
area, the distance from the points chosen in the area to 
Piemenespective atom are calculated in (x',y') coordinates 
and the shortest: distance is chosen. The shortest impact 
parameter is tallied against the atom with which it is 
fets@ctated, in an array. In the case of equal distances 
between two atoms, the atom with which the ion impacts is 
chosen by a uniformly distributed pseudorandom process [17]. 
The final array contains the distribution of all impact 
parameters tallied with each atom in the representative 
area. To obtain the final distribution of impact parameters 
for the crystal face, the individual distributions are 
weighted, added, and the area under the curve is normalized 
fomunity. The number of points used in calculating these 
Mrstriputions are on the order of 100,000 per representative 
area. The final unnormalized distribution of impact parame- 
@ersetor (100) FCC crystals for various angles of rotation 
is shown in Figure 11. Note that all dimensions used in 
@atculatineg these AB aes ons are in terms of Ro» tele 


nearest neighbor distance in the crystal lattice. 
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Figure 11. Unnormalized distribution of impact parameters 
for FCC (100) rotated through various angles 
about <110->. 
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VI. NUMERICAL METHODS AND COMPUTER PROGRAMS 


The theory presented in the preceding chapters has been 
gnicorporated into five FORTRAN IV computer programs which 
are presented in the appendix to this thesis. They appear 
Memeceme Order of their use. Specifically, these programs 
are: 

1. The Herman-Skillman radial-electron density progran, 
Z. The Harrison potential progran, 

fee Distribution of impact parameters progran, 
feeenumerical interpolation program, and 

me ihe Harrison KSE program. 

iomealenlate the coefficient of ion-electron emission 
mr coetficient), these programs are utilized in the fol- 
Mowing manner. The ion-atom-radial-electron densities are 
Galculated using the Herman-Skillman computer program [16]. 
These densities are put into the Harrison potential program 
which calculates the interatomic potential, V(r), used in 
the calculation of the inelastic energy transfer, ETRE" 
The distribution of impact parameters program calculates 
the distribution of impact parameters for each crystal face 
Mig@er investigation.: Since this is a statistical process, 
the distributions recovered from this program are "smoothed" 
uSing a first and tenth degree Legendre orthogonal least 
squares interpolating polynomial, which is generated by 


micmnumerieal anterpolation program [18]. Finally the 
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smoothed distribution of impact parameters and potential 
fimetion are inserted into the Harrison KSE program which 
has been modified to include all the changes to the original 
m@eonry presented in the "Modifications" chapter of this 
mnesis. This program performs the actual calculation of 

meen obs coefficient using the HCM model. The input data 
format to each of these programs is shown in Figures 12 
mereough 16. The definitions of the input variables used 

mm the programs are presented in Appendix A to this thesis. 
mie results presented in the following chapter were run on 


the Naval Postgraduate School IBM/360 Computer System. 
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Data Input for the 
eiectron Density Program 


NNLZCI),WWNLCI), EECI) 


Z,NCCRES ,NVALES , ICN 


RU2(1) 


KEY , TCL, THRISH, MESH, IPRATT , MAXIT ,NOCOPY, KUT 


HEADER CARD 





mepice 12. Yeck setup for rufming the Electron Density 
Program, 
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Patainput tor the 
Potential Program 


ELECTRON DENSITY PROGRAM - OUTPUT ATOM B 


ELECTRON DENSITY PROGRAM - OUTPUT ATOM A 


HEADER CARD 


HNC KSTEPZAZ, AL, BZ,B1,DM,RM,RIEST ,DD,FRR,DRS, IPN 


Figure 13. Deck setup for running the Potential Progran. 
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ICODE , XR, YR,ALPHA, RMNAX 







A6(I,J) 





SCALEF ,RNA,RCL 






ICODE,A,B, INUMBR 







Data Input for the 
Impact Parameters 
Progr ear 


mere 14> Deck setup for runmaing the Distribution of 
Iiodet Panamerers Program. 
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OUTPUT FROM THE DISTRIBUTION OF 
IMPACT PARAMETERS PROGRAM 


HEADER CARD WITH 48 CHARACTERS OF 
INF ORMAT ION 


HEADER CARD WITH 48 CHARACTERS OF INFORMATION 


NPOINT, IF LAG 


Data Input for the 
Interpolation 
eOdnanm 





Hugire 15. Peck setup for rugning the Interpolation Program. 
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Data Input for the 
=e ir @craln 


ZB,ZT,BMAS,TMAS , EOS , KEOM,AO,LSKIP, 
NM,RES ,DRE, REM, IH2 


PSEC) Far(i) 


JF IT 


PACi Le= 1 0) 


EACii ee L = 11s) 


OUTPUT FROM POTENTIAL PROGRAM 


HEADER CARD 


Gu00 0110), (111) 
LCRYST , LTPE, BMAX(1),PUD(1),AROC(I),LMC(1L),KM,FACL 





Figure 16. Deck setup:for running the KSE Program. 
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Vie sesULTS 


Three separate cases of secondary electron emission 
were investigated using the modified HCM model. These were: 
ie Angularly dependent secondary electron emission 
PiometicwlWUG) stace of a copper single crystal 
metated about a <110> axis; 
weeeeoecondary electron emission from Ar, Ne, Kr and 
Mertons mormally incident on the (100), (110) 
eudeqililjyetaces of Cu, Ag, Mo, and Al single 
crystals; and 
eee oecondary electron emission from Ar and Ne ions 
feviialty dncident on the (100), (110) and (111) 
faces of KCl. 
itiecne HCM theoretical model, the primary ion is assumed 
to be neutralized before it strikes the target surface. The 
target atoms in the metal are assumed to be ion cores in a 
Sea of valence electrons. For these reasons, the actual 
results presented here are obtained by modeling the binary 
SelltSion process as a collision between a neutral, inert 
gas atom and an ion core of the target atom. In the cases 
SeecoOpper and silver, the ion cores are fairly well defined, 
eee both appear in column I.B. of the periodic table. 
However, in the cases of molybdenum and aluminum which 
Bereatein the Vl.B: and IJ].A. columns of the periodic 
faible, respectively, there are several possible ion core 
Sent’ curations that could be used in the HCM model. For 
this reason, various degrees of ionization of the molybdenum 


and aluminum target atoms were used in computing the SEE 


coefficient for these target atoms by way of the HCM model. 
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mmeene cases Of copper and silver, results indicate 

miaee the best ion core configuration for both target 
Mererials 1s the singly ionized state. Note that all graph 
titles refer to the HCM model used, whereas the figure 


Gaptions indicate the physical system actually investigated. 


fee ANGULARLY DEPENDENT RESULTS 

The results obtained from the bombardment of Cu single 
crystals by Ar ions for various angles of rotation about 
the <110> crystallographic axis are shown in Figure 17. 
ieeets interesting that the HCM model did not predict the 
meunect behavior of the SEE coefficient for monocrystalline, 
mieubarly rotated targets, but did accurately reproduce the 
mescults obtained for polycrystalline targets [19]. This 
memento indicate that the single collision model is not an 
@ecurate one for arbitrary rotations about the crystallo- 
Paapnic axis, but needs to be modified to account for a more 


eanolicated interaction with the crystal lattice. 


Peete tAL TARGET RESULTS 

An attempt was made Mmctitscmenesss LO Verify those 
Saleulations performed by HCM in their original paper. Not 
Priomsurpriscingly, the results obtained here differ appre- 
erably from those obtained by HCM*using their original 
model. Each system is discussed in detail below. 

1. Copper Bombarded by Ar, Ne, Kr and Xe Ions 

Time i ceorcoinal paper, the Ar - Cu system was 


treated in detail. Using their original model, they 
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° hReCmecppet Crystal. Direction of rotation was . 
about <110>., 
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obtained the results shown in Figure 7, which show excellent 
Meecement with experiment. After altering their model, 
wemne the modifications presented in Chapter IV, the results 
do not agree so well with experiment as did those obtained 
PyeicM. (See Figures 18 through 21.) 

Mimateattempt sta explain this discrepancy, the 
difference between the experimentally observed SEE coefficient 
faoethat obtained here was plotted as a function of energy. 
Meeerecvacally, if the results obtained in this present in- 
Mestigation are to be physically meaningful, this difference 
Must be due to the potential ejection mechanism operating at 
iemier energies. Experimentally, Medvel et al. [1] have 
@eeertained that for Ar bombardment of polycrystalline Mo 
targets, the potential ejection mechanism is still operating 
ee meroics Up to 2.5 Kev, and the PSE contribution is still 
Memdineg to incréase with increasing energy. If this is in 
Bact the case, the results obtained here have physical 
S@emiticance. (See Figures 22 through 25.) 

reais interesting to note that the PSE contribution 
meomeadll three faces investigated exhibit similar behavior 
mm that they all peak at approximately the same energy, and 
Momdetcmdecrease at higher energies. The exact location of 
miempedk, Nowever, Varies With the primary ion type. Further- 
Nene wim all four primary ion cases, the PSE contribution 
Sumves tor the (100) and (110) faces cross at approximately 
Peven tO nine Kev, indicating that the PSE mechanism may 


be somewhat independent of ion type, while showing some 
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dependence on the plane being bombarded. Such a crossing 
Mewavior in the PSE contribution as a function of energy 
mea not be expected simply by analyzing the SEE coefficient 
Sunves, Since these curves do not cross at any energy. If 
fmjes 1s the case, then it would follow that the theoretical 
results obtained for each crystal face should be indepen- 
memely fitted to the experimental data, without regard to 
the other planes being bombarded. In most cases this would 
yield results that are in closer agreement with experiment 
mieyet dO not Climinate the PSE contribution to secondary 
€lectron emission entirely. 

2. oillver Bombarded by Ar, Ne, and Kr Ions 

The silver systems also indicate the possibility 
Semcene PSE mechanism being dependent on the crystal face 
being bombarded as well as the primary ion type. (See 
feoures 26 through 31.) 

In the Ar. - Avesyotenecne Pok contribution to 
Beeendary €lectron emisSion 1S initially largest for the 
(100) face, followed by the (111) and (110) faces, respec- 
miyedy. In all three cases the PSE contribution decreases 
slightly at higher energies. 

In the Kr - Ag system a different ordering of 
itcia) PSE contributions appear. Initially, the (111) 

Pace Contributes the most to potential emission, followed 
Paenee (lO) and (110) faces, ian that order. At higher 
Emergies, the ordering changesy with (100) contributing 


iivemlargest amount to PSE, followed by (111) and (110), 
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Magure 26. Ag bombarded by Ar ions normally incident on the 
CeO 10) and (#11) faces. 
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Mimcetively. Once again, the PSE contribution tends to 
decrease at higher energies. 

In the Ne. - Ag system an entirely different picture 
mrtses. For all energies, the (110) face contributes the 
mist tO PSE, followed by the (111) and (100) faces, 
respectively. However, even at ten Kev, the potential 
Sasson for the (110) face is still rising whereas, over 
mimentire energy range, the PSE contribution from the 
olily and (100) face remains fairly constant. These results 
would seem to indicate, once again, that not only is poten- 
fal emission dependent on primary ion type, but also on 
emyvstal face. 

In the Ar - hems ysStem. tie PSE contribution from 
Pieei lil) face goes to a negative value at seven Kev. The 
meason for this behavior can be found by analyzing the 
Ar - Ag SEE coefficient versus energy curves shown in 
mmere 20.- ihe SEE coefficient experimental data for the 
fll) face tee downward at seven Kev. This behavior is 
not demonstrated on the other faces and would tend to indi- 
Cate that some experimental variation may be present in 
Piese experimental data. On the other hand, the interaction 
being modeled may be sufficiently complicated that the 
potential contribution cannot be separated simply by scaling 
the results. No other system investigated displays similar 


mesults. 
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3. Molybdenum Bombarded by Ar, Ne, Kr and Xe Ions 

inet icures.s2 throuch 43, Mo” initially seems to 
mesthe ion core model that best represents secondary electron 
emission from various target faces. This conclusion iota 
true for the Ne. - Mo and Kr’ - Mo systems. However, for 
the Xe - Mo Seo eiweticebpDest Von core model is not Mo” but 
Mo". This would seem to indicate that the primary ion 
type plays an important role in secondary electron emission. 

4. Aluminum Bombarded by Ar, Kr, and Xe Ions 

In the case of aluminum bombarded by Ar’ (Figures 44 
emrough 52) the Al’ and Al” ion core models for the target 
atoms yielded results that were almost identical. This 
Byctem 1S contrasted with the Carl system which shows a 
Memeinct preference for Al. as the ion core model, with the 
Al and Al. Vicl@dangesiialar and uniformly poor results. 
Similarly, the Mele Al system shows the same preference 
for the Al. ion core model and demonstrates a poorer fit 
when the a and Al’. models are used. When all cases are 
considered, it would seem that secondary electron emission 
iS sensitive to the primary ion type, as has been demonstrated 


by previous results. 


fee ALKALI-HALIDE RESULTS 

Experimental results for the bombardment of alkali- 
halide single crystals by inert gas ions have been obtained 
meomeexperiments by Baboux and Pedrix [21, 22]. Two cases 


are considered in this present™investigation. These include 
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Figure 36. Effect of varying target atom ionization in the 
HCM model on secondary electron emission from 
the Ne - Mo system. 


oS 





ame 


(OMAN DOO OC 


Ne"-Mo 
FE] (111) 


SEE COEFFICIENT 
aS 


NR WW 


pee 





12345678910 
ENERGY (Kev) 


mere o/.  bitect Of Varying target atom ionization in the 
HCM model on secondary electron emission from 
the Ne - Mo system. 


86 





AIP SINAC 
S [1 (110) 


Se COE enema 





12345678910 
ENERGY( Kev) 


mugtren so.  Ettect Of varying target atom ionization in 
the HCM model on secondary electron emission 
from the Kr - Mo system. 


87 





Ss 


Kr-Mo 
r) (111) 


NO 
Ol] 


SEE COEFFICIENT 
Ol 


Os 





emo a > 7 oC 910 
ENERGY (Kev) 


Pecure oo. suateet Of Varyingetarget atom ionization in the 
HCM model on secondary electron emission from 
tiewkr —- Mo system. 


8§ 





ss 


Kr--Mo 
A (100) 


re 
Ol 


SEE COEFFICIENT 
Ol 


O 
Ol 





eee & 7 8 910 
ENERGY (Kev) 


Figure 40. Effect of varying target atom ionization in 
the HCM model on secondary electron emission 
from th. Kr - Mo system. 
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Figure 43. Effect of varying target atom ionization in the 
HCM model on secondary electron emission from 
the Xe - Mo system. 
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Pegvure 44, Etfect of varying target atom ionization in the 
HCM model on secondary electron emission from 
the Ar - Al system. 
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Pigure 45. Effect of varying target atom ionization in the 
HCM model on secondary electron emission from the 
Ar - Al system. 
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meoure 47, Effect of varying target atom ionization in the 
HCM model on secondary electron emission from 
item  - Al system. 
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Eicect sO. Varying target atom ionization in the 
HCM model on secondary electron emission from 
the Kr - Al system. 
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mesure 49. Eitfect of varying target atom ionization in the 
HCM model on secondary electron emission from 
the Kr - Al system. 


98 





SEE COE ean 





[oeo4 5678 910 
ENERGY(Kev) 
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HCM model on secondary electron emission from 
the Xe - Al system. 
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the bombardment of potassium chloride by argon and neon 
mensenormally incident on the (100), (110) and (111) faces 
@rethe single crystal. These two cases are discussed in 


detail below. 
1. Potassium-Chloride Bombarded by Ar Ions 


iiitemcaleulations of the SEE coefficient for alkali- 
mrmeade single crystals present a theoretical problem 
that was not addressed when the metal target results 
werespresented. Specifically, the alkali-halide crystals 
jmeaquestion have lattice structures comprised of two dif- 
ferent atoms arranged in a similar fashion, but displaced 
Mmopace by a distance (1/2, 1/2, 1/2). Because of this, 
jae Standard FCC distribution of impact parameters does 
not work in the alkali-halide case and must be replaced 
yea distribution representative of the normal impact 
Memento With the target surface. The representative 
@reas and distribution of impact parameters are chosen 
mea Manner Similar to that for the metal case, but are 
meronted to include preferential impact with the larger 
mite two 10m Cores in the lattice. The representative 
feeds and distribution of impact parameters are shown in 


mesures 53 through 57. 


The results for KCl bombarded by Ar are shown in 


Figure 58. Experimentally, it was found that the SEE 
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Me@enticient for the (100) face crossed that of the (110) 
meen the five to six Kev region. This behavior, as 
membeas that for the (111) and (110) faces, is predicted 
Meee de HCM modified model. A PSE contribution is evident but 
mas not been investigated. 

2. Potassium-Chloride Bombarded by Ne Ions 


In the case of KCl bombarded by Ne’, the HCM 


modified model predicts a crossing of the (100) and (110) 
MemraceS OLE coefficients as a function of energy at 
approximately five Kev. Such a behavior is not confirmed 
experimentally. 

In an attempt to explain this anomaly, two theories 
mmose., First, such a behavior can be explained, Te ee ne ak Gar 
meee tends to cleave on the (100) face with all the 
chlorine or potassium atoms up. Second, the crossing be- 
havior might be due to preferential sputtering by the primary 
ion of one of the two types of atoms comprising the (100) 
ates leaving a surface that contains either all potassiums 
@r all chlorines. 

inet he tee ase, the crystal structure of the (100) 
face should be the same for both Ar’ and Ne" bombardment, 
mmroicmas prercrential sputtering does not occur. Hence, 
SEE from this face should be due to the bombardment of 
etier pure chlorine or pure potassium, whichever is left 
ime the “up position. 

In the second case, however, preferential sputtering 


: +, 
©f chlorine atoms by Ar is the case most likely to occur, 
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Eamece both the 10n and target atom are of comparable mass. 
iieenhis Situation, preferential sputtering would affect the 
Ar - KCl system more than it would the Ne - KCl Sye tem: 
Since the Ar. - KCl yer Ss SEB cociricient as a function 
Seemeteney 1S accurately predicted for all surfaces examined, 
memdppeans that Sputtering 1s not directly responsible for 
fjresresults that are obtained from either system, 

aie etfect of an alll-potassium or all-chlorine 
Meme cimation to the KC] (100) surface on the SEE coefficient 
‘as a function of energy and primary ion type is shown graphi- 


emty in Figures 58 and 59. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


iiise thesis has attempted to predict the behavior of the 
meetficient of ion electron emission for various primary ion- 
target combinations using the Harrison, Carlston and Magnuson 
theory of secondary electron emission. Three aspects of 
Secondary Clectron emission were investigated, including 
pmeularly dependent secondary electron emission from the 
f100) face of a copper single crystal rotated about <110>, 
secondary electron emission from Ar, Ne, Kr and Xe ions nor- 
meeeeinerdent on the (100), (110) and {111) faces of Cu, Ag, 
Mo and Al single crystals, and secondary electron emission 
from Ar and Ne ions normally incident on the (100), (110) and 
(pei faces of KCl. 

According to the HCM theory, angularly dependent results 
do not demonstrate the non-monotonic dependence on the angle 
of incidence found experimentally. Even so, the HCM theory 
Moeseaccurately model the dependence of the SEE coefficient 
of angle of incidence for inert gas ions normally incident on 
a polycrystalline metal surface. For this reason, it appears 
that the single collision model is not sophisticated enough 
andsturther modification must be made to it, if secondary 
electron emission from rotated, single-crystal faces is to be 
predicted. 

In the case of metal targets bombarded by inert gas 
Honceeitewas totind that the SEE coefficient as a function of 
energy, predicted by the modified HCM movel, deviated from 
that experimentally obtained. This was contrasted with the 
results obtained using HCM's original model which yielded 
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reasonable results in almost all cases investigated. Since 
Piesmodifications presented here were refinements to the 
Original HCM model, it was hoped the results obtained would 
more accurately represent the secondary electron mechanism. 
For this reason, the difference between the SEE coefficient 
Obtained here and the experimental result was identified as 
the potential secondary electron emission component operating 
at higher energies. Such a hypothesis is substantiated by 
experiment [1]. Analysis showed that the potential component 
mimeeteral rose to a peak and then decreased at higher energies. 
Furthermore, this component of secondary emission seemed to 
demonstrate behavior that could possibly be attributed to 
meailary 10m type as well as the face of the crystal being 
bombarded. The reason for this behavior was not positively 
miemearied, but the PSE process appears to be more complicated 
than was originally thought. 

In the. case of potassium-chloride bombarded by Ar ions, 
PaewnGM model fairly accurately predicted the SEE coefficient 
in the zero to 20 Kev energy range for all faces examined. 

On the other hand, potassium-chloride bombarded by Ne ions 
demonstrated a crossing of the SEE coefficient for the (100) 
mam lO )staces as a function of energy which was not predicted 
experimentally. Two hypotheses were proposed that might pos- 
Bio mex larmm this behavior. The first hypothesis proposes 
that KCl cleaves with either all the potassium atoms or all 

the chlorine atoms up, essentially presenting a uniform 


layer of pure potassium or chlorine to the bombarding ion 
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beam, This hypothesis slightly improved the SEE 
coefficient for the Ne’ KCl (U0 Msuniace,~but adversely 
affected the results obtained for the Ar - KCl system. 
The second hypothesis assumes that the (100) surface of Kel 
me~ieecemonstrate preferential sputtering. However, if 
this were the case, preferential sputtering should have been 
most evident in the Ar - KCl System.) Decause Of the Compara— 
Pile smass of Ar and Cl. Since this system demonstrated the 
correct behavior without invoking the preferential sputter- 
Mmiesassumption, it was surmised that preferential sputtering 
loeiot play an important role in the SEE process from the 
meee 100) surface. 

inmaliiecases tested in the present investigation, the 
Pemeoctticient as a function of energy tended to become a 
Memecr approximation to the experimental results at higher 
emeraies. For this reason, it is recommended that primary 
Home careet atom combinations be run for primary ion energies 
between ten and 20 Kev. If in fact the potential component 
of secondary electron emission is real in this energy range, 
it should reach some constant value as primary ion energy 
is increased. Furthermore, because results obtained using 
the modified HCM model are dependent on the ion core configu- 
Prieto meOr the target atoms, it is recommended that more 
work be done in testing the influence of ion-target atom 
EOre lonization oh secondary electron emission. Typically, 
for secondary electron emission from elements such as Cu 


AndwAceewhich appear in column®] of the periodic table, the 
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monmecore clectron configuration is uniquely determined by 
considerations of the atomic valence, but the effect of the 
PmeectrOnic structure of the ion core upon KSE for Mo and 


Al needs more investigation. 
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APPENDIX A 


COMPUTER PROGRAMS 


iMetiis dppendix, the parameters that appear in the 


Someuter programs are presented and explained. 


Variables 


that appear in the data input as constants are stated as 


Suen: 


A pictorial representation of how these data are 


Peieerea into the programs is shown in Figures 12 - 16. 


A. HERMAN SKILLMAN RADIAL ELECTRON DENSITY PROGRAM 


GARD # 


1 


BDO BDO BDO DY DO DO DO PO 


WN 
' 
> 

% 


DmWWW 


2 (Qe) 


NAME 


FORMAT 


DESCRIPTION 


HEADER CARD - CONTAINS 80 CHARACTERS OF INFORMATION 


KEY 
TOL 
Poke SH 
MESH 
TPRATT 
MAXIT 
NOCOPY 
KUT 


RUZ CL) 


Z 
NCORES 
NVALES 
ION 


NNLZ (1) 
WWNL(1) 


EE(1) 


DESCRIBING THE CASE TO-BE RUN 


(14) SEETReterence ( ), 0 

(F8-.6) A PLOO UOT 

(14) " 521 

ais) ‘ » 1 

(14) 2 5 att 

(14) “ » O 

(14) * » 0 

Cae i » 0 

CIPEaS/ Normalized seli-consistent poten- 

Peon s/ jetral as a function of X. Taken 
directly from Reference 16. This 
imam Staueer for cenerating the 
wavefunctions for the atom under 
IMmVvesSt1 eat lon, 

(4. 0) Atomic number of atom in question 

(14) Noe ot shells of atom in question 

(on) 0 

(4) Degree of ionization of atom in 
question 

(14) il iiedtiantum no. Lor Shells 

(F4.0) No. electrons in shell 

(ES. 4)) irialiwercenvalue tor each shell, 


taken from Reference 16. 
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B. POTENTIAL PROGRAM 


CARD # 


i ee 


a 


NAME FORMAT DESCRIPTION 
INUC (12) 0 
KSTEP (13) 10 
AZ (F3.0) REOMLe NO. Ot the first atom 
following the header card below. 
Al (F2.0) lonmzatvon Of the first atom 
following the header card. 
BZ (E520) Atomic no. of the second atom 
following the header card. 
BI (E220) Jlonization of the second atom 
following the header card. 
DM (E220) 5 
_ RM (C2 0 ) 5 
REEST (2-0) 5 
DD (5). ly) .020 
FRR Cro.) ry 
DRS CES. 1) SOan0 
IPN (14) 1 


HEADER CARD - CONTAINS 30 CHARACTERS OF INFORMATION 
PESCRURING THE “CASE (0) BE RUN 


OUTPUT EROMSELECTRON DENSITY PROGRAM WITH FIRST 
CARD REMOVED CORRESPONDING TO THE ATOM DESCRIBED 
BY AZ AND AI ABOVE 


CULPUT FROM® ELECTRON DENSITY PROGRAM WITH FIRST 
CARD REMOVED CORRESPONDING TO THE ATOM DESCRIBED 
BY BZ AND BI ABOVE 


C. DISTRIBUTION OF IMPACT PARAMETERS PROGRAM 


CARD # 


1 


NAME FORMAT DESCRIPTION 


HEADER CARD - 80 CHARACTERS OF INFORMATION 
DESCRIBING eInE CASE TO) BE RUN 


ICODE CIs) Ie tierepresentative area 15 
a rectangle 
2, if representative area is 
Ae ie lee lee 


A (E10. 4) HeimgibeoLteatepbesemtative area 
i eKOmuUn1 tS. 

B (F10.4) Length of representative area 
Pee Unt es 

INUMBER (13) lem ot atomssin representative 
area. 
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feeb iLolRIBUTION OF IMPACT PARAMETERS PROGRAM (Continued) 


CARD # NAME FORMAT DESCRIPTION 


5 EARGE CT) (12) bogvealearray that indicates if 
any particular atom is larger 
Cianeerhemotier, ise. , Cl in 
KO seS 82.9 no, It yes 
CNeEsdIStribution 1S) skewed to 
iidhedrteuptenrercntial collision 
with the larger atom. 


4 SCALEF (RIL a) 1 g10: 

4 RNA (F10.4) Tonic radius of smaller of two 
atoms. 

4 RCL (F10.4) Ionic radius of larger of two 


atoms--Both RNA and RCL can 
be set to anything if atoms 
are both the same size. 


5-7 AS C1. J) (Slee nearnay that contains thie 
Gey) evereices Of a triangular 
representative area in Ro 


UViis elt inewared 15 nOt tri 
aneurin tiey can be Set to 
anything. 


8-11 NOE) ) ile eit arnayetnat Comtains the (x,y,z) 
Position of atoms in the repre- 
sentative area in Ro units. I 
runs from 1 to INUMBR, J runs 
frome ls co o.. 


eZ Pit (P1024) VAmelesan radians x'-axis makes 
with x-axis. 

a2 XR : CEO 4) eee posit lon soLe(x',y') origin. 

12 YR C024 iy position. of (x!,y") origin. 

2 ALPHA (EGE 4) Seinele sin radians che representa- 
tive area is rotated through. 

eZ RMAX (F10.4) Maximum impact parameter 


expeeted. The mo. 25) is 
WSs Ud Pre pent. 


is 





D. NUMERICAL INTERPOLATION PROGRAM 


CARD # NAME FORMAT DESCRIPTION 


i NPOINT (14) No. points from the distribution 
of impact parameters program to 
be smoothed. 

i IFLAG (ra } iPS indreates writtem output only. 
Mernavucares written. Output), 

card output and graphical 
output desired. 


iz HEADER CARD - 46 CHARACTERS TO BE OUTPUT WITH GRAPHS 
S HEADER CARD - 46 CHARACTERS TO BE OUTPUT WITH GRAPHS 
4 OUI SPROMeDISTRIBUTION OF IMPACT PARAMETERS PRO- 


GRAM WITH ALL "NEW ATOM" CARDS REMOVED 


Eee kKSE PROGRAM 


CARD # NAME FORMAT DESCRIPTION 
i eR Yo T (A4) Crysraletype. FCC!/ BCC" 
It PE (iz) 1 
i BMAX (1) (3F6.4) Maximum impact parameter in Ro 
Win GeeerOr cach face . 
PUD(1) (C5E6 4) Ry fHOoreecdehe face im RO Units. 
ARO(T) (3F6.4) "Nos atoms per representative 
; Dreamin eROmunits, for each face, 
] LM(I) (3135) Totalenumeer Of points in each 
dvstriabucrom, 
1 KM 4) 250 
1 FACL (F7.4)  (2)1/4/4 for Fcc, (3)1/%/4 for 
Bee 
2-A PEP Cr) (2014) Distributton Of Impact parameters 
Proimwinterpolatvon Program. 
B IH1 80A1 SUncitaedeters of intormation. 
C 1S 19A4 fo characters Of intormation from 
fivrsuscard sot voueput trom Potential 
Program. 
C KMR (14) 250 
D APSI (I) ILS Wave functions from the Potential 
program, Atom A. 
D BPSI (TI) OES. 5 Waue functions from the Potential 


pmogram, Atom bE. 
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E. KSE PROGRAM (Continued) 


CARD # 


NAME 


FORMAT 


DESCRIPTION 


ep) 


epeq pas yanaq aq ac i) (86) (Se) (ae, Sm jae, jae) 20 


Vay) 


EA(1) 
IAH 


EA(I) 
Sient 


Poe cl) 


Elr (J) 


ZB 
aM 
BMAS 
TMAS 
EOS 


DEO 
KEOM 


AO 
lk IP 
NM 
RES 
DRE 
REM 
IH2 


e245 


8F8.4 
4A4 


CES. 5) 


(F5.4) 


Potential between atoms A and 
BerLonericmeocentlal Frogram, 


SpPecthescopie 21On1 zation 
potentials, numbers 11 to 18. 
MCOmelanaetersvort accounting 
information. 


spect rosecopie 10Nni zation 
POten tals enumpers 1 to 10. 
No. of values to which results 
ages tants 

YpsE experimental for curve 


being efi k. 
Y¥eexperimental for (lll) face 
Ate eke ve 
ALOMUC nO tor atom A. 
AgomiceNo. tor atom B. 
Atomic weight for atom A. 
Atomic weight for atom B. 
Energy in Kev at which resuits 
aLemeOustarl., 
EVeio ves cen ain Kev. 
Energy in Kev at which results 
age ste end: 
Pe mlecrecOons tant sin angstroms. 
il 
cS 
50 
50 
5000 


Accounting data. 


Smidmeatescsthat the actual cardynumber varies from run to 
mineandenas been replaced by @ letter to indicate the 
order in which the cards are arranged. 
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Misismtnesis 1S a computer simulation of secondary electron 
emission (SEE) from monocrystalline metals and alkali-halides 
Msane a modified version of the Harrison, Carlston and Magnuson 
menele collision theory of SEE. 
Bated: the angular dependence of SEE from Cu bombarded by Ar , 
the dependence of SEE as a function of energy for rare gas 
mise normally incident on the (100), (110) and (111) faces of 
Secal single crystals, and the dependence of SEE as a function 
Seeenergy for Ar and Ne ions normally incident on the (100), 

The theory does not accurately 
Mescribe the angular dependence of SEE for monocrystalline Cu 
mireets, but does accurately predict the modified sec 6 dependence 
Bouma experimentally in polycrystalline studies. For the metal 
targets, the difference between the theoretical kinetic secondary 
memtesion result and the experimental datum is identified as 

The alkali-halide SEE simulation 
agrees reasonably well with experiment. 


mero) and (111) faces of KCl. 


potential secondary emission. 
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